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ABSTRACT 

It  is  evident  that  military  applications  for  21s'  century  will  be  highly  complex,  multivariable  systems. 
Designing  optimal  controllers  for  such  applications  will  require  the  use  of  mathematical  models  which 
describe  the  complexity  of  the  underlying  processes  as  accurately  as  possible.  For  robustness  it  is  essential 
that  these  models  include  the  uncertainties  in  the  estimated  process  dynamics  and  trajectories.  Controllers 
resulting  from  optimal  control  strategies  (like  LQG,  H2,  Hf  will  usually  be  of  very  high  order  and  that  can 
cause  implementation  and  computational  problems.  Some  of  these  problems  can  be  overcome  by  using 
controllers  that  are  of  lower  order  and  restricted  structure.  However,  the  design  and/or  tuning  of  such 
controllers  in  order  to  provide  the  performance  comparable  with  full-order  solutions  is  still  a  very 
contentious  issue.  Added  to  this  is  the  need  to  achieve  robust  properties  and  performance  specifications 
required  by  military  applications.  Very >  few  methods  for  designing  restricted-structure  controllers  exist  that 
allow  the  performance  and  robustness  objectives  to  be  combined  into  one  relatively  simple  optimisation 
problem.  This  lecture  presents  an  LQG/H2-based  method  that  tackles  the  above  mentioned  issues. 
The  LQG/H2-based  criterion  is  minimized  in  such  a  way  that  the  resulting  controller  is  of  the  desired  form 
and  is  causal.  A  simple  analytic  solution  cannot  be  obtained,  however  a  straightforward  direct  optimization 
problem  can  be  established  which  provides  the  desired  solution.  From  this  solution  it  is  also  possible  to  derive 
a  tool  for  assessing  the  performance  of  existing  controllers. 


1.0  INTRODUCTION 

The  design  and  performance  assessment  of  optimal  controllers  (Kwakemaak  and  Sivan,  1972  [1],  Desborough 
and  Harris  1993  [29],  Grimble  and  Uduehi  2001  [26])  is  considered  as  well  as  the  ways  in  which  robustness 
properties  can  be  modified  and  assessed.  The  minimization  of  quadratic  cost  functions  for  stochastic  linear 
systems  in  such  a  way  that  the  robustness  margins  are  improved  is  of  interest.  There  are  of  course  well  known 
guaranteed  robustness  properties  for  systems  with  state  feedback  (Anderson  and  Moore  1971  [2],  Safonov  and 
Athans  1986  [3],  Doyle  1978  [4])  and  in  some  problems  these  properties  can  be  recovered  using  a  loop 
transfer  recovery  approach  ([5]-[10]).  However,  more  direct  methods  of  tuning  robustness  properties  would  be 
valuable  (Doyle  and  Stein,  1981  [11]).  This  is  explored  for  both  full-order  optimal  output  feedback  controllers 
and  for  low  order  (restricted  structure)  optimal  solutions. 

This  paper  also  examines  the  important  issue  of  performance  assessment  of  existing  controllers  with  regard  to 
robustness  and  other  steady  state  properties.  In  recent  years  there  has  been  a  lot  of  research  devoted  to 
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assessing  the  performance  of  existing  control  systems  with  the  objective  of  determining  if  some  degrees  of 
freedom  exist  to  optimise  such  control  systems.  The  idea  is  to  compare  the  performance  of  the  existing  control 
against  some  known  reference  standard  and  produce  an  index  ranging  between  0  and  1 ,  which  indicates  how 
well  the  control  system  is  performing.  Desborough  and  Harris  (1993  [29])  considered  the  assessment  of 
control  loop  performance  for  both  feedback  and  feedforward  control  using  minimum  variance  as  the 
benchmark  cost  measure.  Huang  and  Shah  (1999  [30])  summarised  the  state  of  the  art  in  a  monograph,  mostly 
focussing  on  the  minimum  variance  cost  index  as  the  performance  assessment  measure.  By  reducing  the 
variance  of  the  output  deviation  systems  can  be  operated  as  close  to  their  physical  constraints  as  possible. 
However,  due  to  the  problems  in  implementing  minimum  variance  controllers,  ( McGregor  and  Tidwell  1977 
[31],  wide  bandwidth,  large  noise  amplification)  and  because  the  minimum  variance  controllers  tend  to  be  of 
high  order  relative  to  classical  designs  (such  as  PID),  the  issues  of  performance  assessment  were  extended  to 
include  LQG,  GPC  and  GMV  criteria  for  restricted  structure  reduced  order  controllers  (Grimble  [32],  Ordys 
et  al  [28],  Uduehi  et  al  [27]).  However,  none  of  the  performance  assessment  indices  developed  so  far  tackle 
the  issues  of  robustness. 

Two  aspects  of  the  design  of  full-order  and  restricted  structure  optimal  controllers  are  considered  in  the 
following.  The  first  involves  robustness  improvement  and  the  second  is  concerned  with  the  noise  rejection 
properties.  An  LQG  criterion  is  to  be  optimized  but  the  usual  robustness  and  noise  rejection  properties  will  not 
hold  if  the  controller  structure  is  limited  to  say  a  PID  or  a  low  order  lead-lag  form.  Grimble  (1999  [12], 
2000  [13])  introduced  a  polynomial  systems  approach  to  restricted  structure  optimal  control  design  and  this  is 
the  philosophy  followed  here.  However,  the  design  of  such  controllers  previously  focussed  on  performance 
issues  and  the  robustness/noise-rejection  aspects  were  not  considered  in  any  detail.  The  strategy  for  robustness 
improvement  is  to  add  a  fictitious  signal  and  a  sensitivity  costing  term  in  the  criterion.  This  enables  the 
penalty  on  sensitivity  to  be  directed  at  modifying  the  robustness  properties.  The  normal  LQG  cost-index  does 
of  course  include  error  and  control  signals  that  depend  upon  the  sensitivity  functions.  However,  this  does  not 
enable  these  sensitivity  terms  to  be  costed  in  a  particular  way.  The  proposed  robustness  weighting  term  gives 
free  choice  of  the  weighting  function  and  enables  the  H2  norm  of  a  weighted  sensitivity  function  to  be 
minimized.  This  relates  to  the  definition  of  a  so-called  Dual  Criterion  (Grimble,  1986  [14]),  but  the  results 
here  are  focussed  on  the  design  issues  and  they  use  what  might  be  termed  a  Kucera  polynomial  systems 
approach,  (1980  [15]).  The  impact  of  a  coloured  measurement  noise  model  on  the  controller  design  and  on  the 
robustness  and  performance  properties  is  examined.  The  frequency-domain  polynomial  systems  approach  is 
particularly  helpful  when  determining  the  frequency  response  behaviour  to  noise  and  disturbance  signals 
(Grimble,  1994  [16]).  The  questions  to  be  answered  are  how  the  sensitivity/robustness  weighting  and 
the  measurement  noise  models  can  be  selected  to  optimize  the  robustness  and  noise  rejection  properties. 
The  impact  of  such  a  weighting  on  both  the  full-order  optimal  and  on  the  restricted  structure  controllers  will 
be  explored. 

2.0  SYSTEM  MODEL 

The  system  shown  in  Fig.  1  is  assumed  to  be  linear,  continuous-time  and  single-input,  single-output. 
The  external  white  noise  sources  drive  colouring  filters  which  represent  the  reference  Wr(s),  measurement 
noise  Wn(s),  robustness  modification  Wp(s)  and  disturbance  Wd(s)  subsystems.  The  robustness  model  does  not 
exist  physically  but  is  introduced  for  design  modification.  The  system  equations  become: 


Input  disturbance: 

d(s)  =  Wd(s)q(s) 

(2.1) 

Robustness  signal: 

p(s)  =  Wp(s)r/(s) 

(2.2) 
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Output: 

y(s)  =  d(s)  +  p(s)  +  W(s)u(s) 

(2.3) 

Reference: 

r(s)  =  Wr(s)C>(s) 

(2.4) 

Tracking  error: 

e(s)  =r(s)  -  y(s) 

(2.5) 

Observations: 

z(s)  =  y(s)  +  n(s ) 

(2.6) 

Measurement  noise: 

n(s)  =  W„(s)a)(s) 

(2.7) 

Control  signal: 

u(s)  =  C0(s)(r(s)~z(s)) 

(2.8) 

The  system  transfer  functions  are  all  assumed  to  be  functions  of  the  Laplace  transform  complex  number  in  the 
complex  frequency  domain.  For  notational  simplicity  the  arguments  in  W(s),  and  the  other  models  are  omitted. 


Disturbance  Noise 


Observations 

z 


Figure  1:  Single  Degree  of  Freedom  Unity  Feedback  Control  System  with 
Measurement  Noise,  Reference  and  Disturbance  Models. 


Assumptions 

1.  The  white  noise  sources,  co,  ij  and  f  are  zero-mean  and  mutually  statistically  independent. 
The  intensities  of  these  signals  are  without  loss  of  generality  taken  to  be  of  value  unity. 

2.  The  system  W  is  assumed  free  of  unstable  hidden  modes  and  the  reference  Wr ,  noise  Wn  ,  robustness  Wp 
and  disturbance  Wd  subsystems  are  asymptotically  stable. 

The  following  expressions  may  easily  be  derived  for  the  output,  error,  observations,  controller  input,  control 
and  sensitivity  costing  signals: 

Output:  y(s)  =  WC0(\  +  WC0y\r(s)-n(s))  +  (l  +  WC0y\d(s)  +  p(s))  (2.9) 

Error:  e(s)  =  r(s)  -  y(s)  =  (1  +  WC0  )_1  (r(s)  -d(s)~  p(s))  +  WC0  (1  +  WC0 )_1  n(s)  (2.10) 
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Observations: 

z(s)  =  WC0(l  +  WC0)~lr(s)  +  (1  +  WCy~\d(s)  +  n(s )  +  p(s )) 

(2.11) 

Controller  input: 

e0  (s)  =  r(s)  -  z(s)  =  (1  +  WC0  )_1  (r(s)  -d(s)-n(s)~  p(s )) 

(2.12) 

Control  signal: 

u(s)  =  (1  +  WCoyl  C0 (r(s)  —  d(s)  —  n(s)  -  p(s)) 

(2.13) 

Robustness  signal: 

h(s)  =  (l  +  WC0y1p(s) 

(2.14) 

These  equations  include  the  following  sensitivity  operators: 

Sensitivity:  S  =  (1  +  WC0 y1  (2.15) 

Complementary  sensitivity:  T  -  1-S-  WC0S  (216) 

Control  sensitivity:  M  =  C0S  =  C0  (1  +  WC0  f1  (2.17) 

The  system  shown  in  Fig.  1  may  be  represented  in  polynomial  form  ( Kailath  1980  [22]),  where  the  system 
transfer  functions  are  written  as: 


System: 

w=a~!b 

(2.18) 

Reference  generator: 

1 

II 

(2.19) 

Input  disturbance: 

Wd  =  A~’cd 

(2.20) 

Measurement  noise: 

Wn=A~lCn 

(2-21) 

Robustness  signal  model: 

Wp=A~1Cp 

(2.22) 

There  is  no  loss  of  generality  in  assuming  these  models  have  a  common  denominator  A  polynomial. 
The  various  polynomials  are  not  necessarily  coprime  but  the  system  transfer  function  is  assumed  to  be  free  of 

unstable  hidden  modes.  The  coprime  representation  of  the  system  is  denoted  by  A0lB0,  where  B=BqUq  and 

A=AoUo . 

The  spectrum  of  the  signal  r(s)-d(s)-n(s)-p(s)  in  equations  (2.12)  and  (2.13)  is  denoted  by  0  ^  (s)  and  a 
generalised  spectral-factor  17  may  be  defined  from  this  spectrum,  using: 

YfY*f  =  0#  =  O,,  +  <$>dd  +  0„„  +  ®pp  (2.23) 

In  polynomial  form  Yj-  =  A  1  Dj  .  The  disturbance  model  is  assumed  to  be  such  that  Df  is  strictly  Hurwitz  and 
satisfies: 

DjD)  =  EE  +  CdC*d  +  C„d  +  CpC;  (2.24) 

The  role  of  the  robustness  model  Wp  may  now  be  explained  since  it  is  one  of  the  components  in  the  combined 
signal  spectrum  0ff  defined  in  (2.23).  Clearly,  if  0pp  is  dominant  then  the  spectrum  0ff  — » 0pp . 
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This  implies  that  <I>  — »  S®ppS*  and  this  spectrum  will  provide  the  robustness  term  (suitably  weighted) 

introduced  in  the  cost  function,  in  the  next  section. 


3.0  LQG  CRITERION  AND  RESTRICTED  STRUCTURE  CONTROL  PROBLEM 


The  LQG  cost-function  to  be  minimized  ( Yoiila  et  al  1976  [17])  is  defined  as: 

J  =  {Qcmee(s)  +  Rc(s) 0,»  +  Pc(.s')Oeoeo (S)}ds 

J  D 


(3-D 


where  Qc,  Rc,  Pc  represent  dynamic  weighting  elements,  acting  on  the  spectra  of  the  error  e(t),  feedback 
control  u(t)  and  controller  input  e0(t)  signals.  The  Rc  weighting  term  is  assumed  to  be  positive  definite  and 
Qc  ,  Pc.  are  assumed  to  be  positive-semidefmite  on  the  D  contour  of  the  s-plane.  The  robustness  weighting 
term  can  be  motivated  as  in  the  dual  criterion  results  of  Grimble  (1986  [14]).  The  error,  control  and  robustness 
weightings  can  be  written  in  polynomial  form  as: 


Qc  = 


AA 


M 

AA 


R,=- 


R, 


B,n 


44  A  A 


and 


P  = 


BpBp 

A\Aq 


(3.2) 


where  Aq  is  a  Hurwitz  polynomial  and  Ar  is  a  strictly  Hurwitz  polynomial  ( Grimble  and  Johnson  1988  [19]). 
The  problem  will  be  to  minimize  the  above  criterion  with  the  controller  chosen  to  have  a  specified  structure 
[13].  The  Pc  weighting  term  represents  a  robustness  weighting  to  be  explained  later.  For  later  use  let  the 
weightings  be  rewritten  using  the  common  denominator,  so  that: 


Qc  ~  Qc  +  Pc  ~ 


O  +P 

AA 


AqAq 


or  Qc  =  Qa,  KAqA  A,. A,.)  =  Qcn  /(AWAW) 


where 


Qcn  ~  QaA-A*  and  Aw  =  AqAr  .  Similarly  for  the  control  weighting: 

A=-yA  =  An  !{AqA*qArA*r)  =  Ren  KAX) 

A  A 

where  Rcn  =  RcnAqA*q  . 

Theorem  3.1:  Restricted  Structure  Single  Degree  of  Freedom  LQG  Control  Problem 

Consider  the  LQG  error  and  control  weighted  criterion  defined  in  (3.1),  and  the  system  introduced  in  §2. 
The  conditions  that  determine  the  LQG  controller  of  restricted  structure  are  derived  below.  The  derivation 
includes  the  robustness/sensitivity  costing  and  coloured  measurement  noise  model.  These  terms  were  not 
considered  in  the  previous  polynomial  approaches  to  the  restricted  structure  control  design  problem.  The  cost- 
function  to  be  minimised  was  defined  in  as, 

A  —  ~  {Qc^ee+RcQuu+Pc®^ 

2*Jd 
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Noting  the  independence  of  the  noise  sources  and  recalling  (2.10)  to  (2.13)  and  the  definitions  for  sensitivity 
in  equations  (2.15)  to  (2.17)  obtain  by  substituting  in  (3.1): 


J  = 


2  ft  j  d 


f  Qc(l-  WM)(X>ff(\  -  W*M*)  +  PC(  1  -  WM)$>  ff(\  -  M*W*) 
\+RcM®  ffM* -Qc(\-WM)®nn-®nn(\-M*W*)Qc+Qc®nn 


j  =  _J_r  \{W*QCW  +  RC)M O  ffM*  -M*W\QC<$>  ff  -  Qc<S>nn)  1 
2*jl  \-(Qc®ff-Qc<S>nn)WM  +  Qc®ff-Qc®nn  J  5 


(3.3) 


where  Qc  -  Qc  +  Pc  and  this  may  be  written  in  the  alternative  forms  Q,  —  ^r":.  =  "  + or  using  a 


AqAq  AqAq 


common  denominator  for  the  spectral  factor: 

Qc  -  Qcn  l(AqA*qArA*r  )  -  Qm  /(AX) 

where  Qcn  =  QcnArA* .  Then  define, 

®fi=W(QcQff-Qc®m) 


(3.4) 


(3.5) 


and  0^=0  +0jj+0  +®  . 

jj  rr  da  nn  pp 

The  generalised  spectral  factors  Yc  and  Y,  due  to  Shaked  (1976  [21])  may  now  be  defined,  using: 

Y*YC  =  W*QCW  +  Rc 

YfY*  =  =  (t>rr  +  ®dd  +  <X\n  +  <$>pp 


Completing  the  squares  in  equation  (3.3)  obtain: 


where 


J  =  {(YtMYf-^XYcMYf-^)*  +  O0}ds 

2njD  vv  vv 


n-  y, 


YcYf 


d>0  =  Qc®ff-Qc<*>nn- 


®fp*jp 

%  H= 

YcYJfYf 


(3.6) 

(3.7) 


(3.8) 


(3.9) 


Substituting  in  the  spectral-factor  expressions  (3.6)  and  (3.7),  using  the  polynomial  system  models  in 
equations  (2.18)  to  (2.22),  obtain 


YfY}  =  (EE*  +  CdC*  +  C„C;;  +  CpC*p)/(AA*)  (3.10) 

Y}YC  =(B*QcnB  +  A*RcnA)/(A*A*wAwA)  (3.11) 
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where 

Qcn=(Qcn  +  Pcn)ArAl  and  Rcn  =  RcnAqA*q  (3.12) 

To  obtain  the  polynomial  spectral-factors  define  the  filter  Dj  and  control  Dc  spectral  factors  from  (3.6)  and 
(3.7)  respectively,  as: 

DjDf  =  EE  +  CdC*d  +  CnC*n  +  CpC*p  (3.13) 

D*DC  =  B*QcnB  +  A*  Rcn  A  (3.14) 

Recalling  that  Aq  and  A r  are  normally  chosen  to  be  coprime,  the  generalized  spectral  factors  may  be  written  in 
the  form: 

Yf-A~!Df  and  YC  =  A~1DC  where  AC=AAW  and  Aw=AqAr 

The  various  terms  in  the  criterion  (3.1)  may  now  be  simplified  by  substituting  from  the  polynomial  system 
models  in  §2  and  the  spectral  factor  results  given  above. 


*>Jp=W  (Qc®ff-Qc®„n)  =  —( 


B*  QcPfDf  -QcnCnCn 


ff 


A 


*  * 

AqAqAA 


) 


O fp/(YcYf)  =  B  ( 


QcnDfDf 

AA„ 


QcnC  nPn_y  A- 


AA„ 


DcDf 


(3.15) 


The  measurement  noise  subsystem  must  be  asymptotically  stable  and  after  cancellation  of  common  terms  may 
be  written  as:  A~lCn  =  A~qCii0  ,  where  A  -  An()  A() .  The  following  diophantine  equations  must  be  introduced: 

Feedback  diophantine  equations: 

Calculate  ( G0 ,  H0,  F0)  with  F0  of  minimum  degree: 

D*cG0+F0AAq=B*QcnA;Df 
DCH0  -F0BAr  =  A  RcnAqDf 

Implied  equation: 

Multiplying  (3.16)  by  BAr  and  (3.17)  by  AAq  and  adding  the  equations  obtain: 

D*c(G0BAr  +  H0AAq )  =  (B*QcnB  +  ARcnA)Df 

•  •  •  ^  • 
and  after  division  by  Dc  obtain: 

G0BAr  +  H0AAq  =  DcDf  (3.18) 


(3.16) 

(3.17) 
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Measurement  noise  equation: 


DcDfX0  +  Y0An0A  =B  QcnCn0CnAr 


(3.19) 


The  terms  in  the  cost-optimization  problem  may  now  be  considered,  substituting  from  the  above  polynomial 
equations.  Substituting  from  the  diophantine  equations  (3.16)  and  (3.17): 


O 


fp 


Y  Y 

V  V 


=  C 


Jo 


+4)-c 

aa^  Dy 


X „ 


AnOAg 


DcDf 


(3.20) 


Considering  now  the  term  Yc  MYj  ,  writing  C0  =  C01dC0n ,  obtain  : 


DcDfC0n 


Y  MY f  =- 

f  AAJAC0d  +  BC0n) 

The  first  squared  term  in  (3.8),  using  (3.21)  now  becomes: 

v  .  ,v  ® fp  [DcDfC0n  ~(G0 -X0Ad)Ar(AC0d  +5Co„)]  ( Y0-F0D y) 

YcMYf~TpXX^ - - - +  ■ 


AAw(AQ)d  +  BC0n) 


* 

DcDf 


(3.21) 


(3.22) 


Substituting  from  the  implied  diophantine  equation: 

v  .  .v  ®  fp  [(HoAnOAq  +X0A,.B)C0n  ~(G0A  0 -X0A)A  C0d]  (Y0-F0Df) 

YcMYf-7FYX  = - - - +  ' 


AnOAw(AQ)d  +BC0lJ) 


% 

DcDf 


(3.23) 


This  cost  term  expression  may  be  written  in  the  form: 

YcMYf-—^  =  Tl++Tl-  (3.24) 

YCYf 

where  the  term  within  the  square  brackets  in  (3.23)  denoted  by  7)'  .  This  term  is  stable,  since  Aw  is  Hurwitz 
and  the  closed-loop  characteristic  polynomial  pc  -  (AC0d  +  BC0n)  is  required  to  be  strictly  Hurwitz  for 
Jmin  <  oo .  The  final  term  in  (3.23)  is  strictly  unstable  since  Dc  is  strictly  non-Hurwitz. 


4.0  COST  FUNCTION  MINIMIZATION  AND  PARAMETRIC  OPTIMIZATION 
PROBLEM 

Given  the  simplification  of  terms  in  the  cost-function  presented  above  the  cost  minimization  procedure  may 
be  followed  ( Grimble  and  Johnson,  1988  [19]).  Note  that  the  cost-  function  (3.8)  may  be  written,  using  (3.24) 
as: 

J  =  {(Y1  +  TJ )(Tq  +  Tf  )*  +  <!>0}ds  (4.1) 
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From  the  Residue  theorem  the  integrals  of  the  cross-terms  7}+7}  ,  7}  7}'  can  be  shown  to  be  zero.  This  result 
follows  because  Q)  Tj  Tj  ds  =  —  d)  7,  7,  <7.s  but  the  term  T} '  is  analytic  for  all  s  in  the  left  half  plane  so 

that  the  sum  of  the  residues  obtained  in  calculating  §  Zj  Tx  ds  is  zero. 

Note  that  this  result  still  applies  if  the  function  Tx  7]+  contains  poles  on  the  jco  axis,  since  they  can  be  avoided 
by  the  D  contour,  using  small  semi-circular  detours  in  the  left-half  plane.  These  semi-circles  are  centred  on 
these  poles  and  do  not  contribute  in  the  limiting  case  as  the  radius  tends  to  zero.  Also  observe  that  the  term 

containing  TfT*  could  lead  to  an  infinite  cost  should  such  terms  be  present.  Flowever,  these  may  not  be 
present,  since  the  optimal  control  may  be  chosen  so  that  they  cancel.  The  practical  case  when  this  arises  is 
when  the  error  weighting  includes  an  integrator  A  (s)  =  s.  When  the  controller  denominator  C0d(s)  includes 
integral  action  the  A  polynomial  cancels  throughout  the  term.  The  consequence  is  that  the  criterion  can  have  a 
finite  minimum,  even  though  certain  cost  function  terms  include  j  axis  poles.  The  cost-function  therefore 
simplifies  as: 

J  =  ^—§  {(W  +  TfTf*)  +  %}ds  (4.2) 

2*  J  D 

Since  the  terms  Tj  and  are  independent  of  the  controller,  the  criterion  J  is  minimised  when  the  first  term 
involving  7)'  is  minimized.  Flowever,  if  the  feedback  controller  C0  has  a  restricted  structure  then  it  is  unlikely 
that  Tj+  can  be  set  to  zero.  It  follows  that  to  minimize  the  cost-function  the  first  term  in  (4.2)  should  be 
minimised,  through  the  choice  of  C0,  namely: 

Jo=^~HTi+Ti+*}ds  (4-3) 

For  a  finite  solution  to  this  cost  minimization  problem  to  exist  the  7}+  term  must  be  asymptotically  stable. 
Inspection  of  this  term: 

j,+  _  [(H0An0Aq  +  X0ArB)C0n  -  (G0An0  -  X0A)ArC0d  ]  ^ 

Ai0^w(dC0d  +BC0n ) 

reveals  that  all  terms  are  asymptotically  stable  but  the  weighting  Aq  could  include  a  j  axis  zero  (A  is  only 

assumed  to  be  Hurwitz).  However,  it  is  assumed  that  although  the  structure  of  the  controller  C0  =  C0„Coj  is 
limited,  C0d  will  have  zeros  at  the  j  axis  zeros  of  the  chosen  weighting  Aq.  Thus,  such  a  zero  will  cancel  and 

under  the  given  assumptions  7}+  is  asymptotically  stable. 

Then  the  LQG  controller  of  restricted  structure  may  be  calculated  from  a  simple  direct  optimization  problem. 
First  compute  the  filtering  and  control  spectral  factors  Df  and  Dc  (strictly  Hurwitz  due  to  the  system 
description)  using: 

DfD}  =  EE*  +  CdC*d  +  CnC*n  +  CpC;  (4.5) 

D*DC  =  B*QcnB  +  A*RcnA  (4.6) 
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The  following  regulating  diophantine  equations  must  then  be  solved  for  (G0,  H0,  F0),  with  F0  of  minimum 
degree: 

D*cG0+F0AAq=B*QcnA*rDf  (4.7) 

D*H0  -F0BAr  =  A*RcnA*Df  (4.8) 

and  the  following  measurement  noise  diophantine  equation  must  be  solved  for  (X0,  Y0),  with  Y0  of  smallest 
degree: 

D>*fX o  +  Y0An0Aq  =  B*QcnCn0C;X  (4.9) 


The  optimal  controller  C0  =  C0nC0d  must  then  be  found  to  minimize  the  following  component  in  the 
cost-function  term: 


where 


J0=^\[  iK (joj)Tt+ (-joj)}dco 


(4.10) 


[(H0An0Aq  +  X0ArB)C0n  ( G0An0  X0A)ArC0d] 
4«o  Ar(A  Cod  +  BC o,,) 


If  the  controller  has  a  specified  limited  structure  the  minimum  of  the  cost  term  J0mjn  will  be  non-zero.  For  an 
unconstrained  solution  the  minimum  is  achieved  when  7)'  =  0  and  the  minimum  of  J0  (denoted  J0min)  is  zero. 
It  has  been  explained  in  Theorem  3.1,  that  the  computation  of  the  optimal  feedback  controller  Cq  reduces  to 
minimization  of  the  term  J0. 

It  is  clear  from  (4.4)  that  Tf  can  be  written  in  the  form: 

V  HCon^-Co^KC^  +  CM  (4. 1 1) 

where  C0  =  C0n  /  Cod  has  a  specified  structure  which  can  be  as  expressed  below: 


Reduced  order: 


„  ,  .  cn0 +cnls+---+cnpSP 

C0(s)  = - r 

cd0  +cdls+---+cdvs 


where  v  >  p  is  less  than  the  order  of  the  system  (plus  weightings) 

Lead  lag: 

q  ssj  _  ( cn0  +  cnls)( cn2  +  cn 3s) 
(cd0+cdls)(cd2+cd3s) 
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PID: 


Co  ( s)  —  kg  +  kj  /  s  +  k2s 

The  assumption  must  be  made  that  a  stabilising  control  law  exists  for  the  assumed  controller  structure. 
Note  that  the  controller  structure  should  be  consistent  with  the  choice  of  error  weighting,  in  the  sense  that, 
if  Aq  includes  a  j  axis  zero,  then  the  controller  denominator  Cq^(s)  should  also  include  such  a  zero. 
The  solution  of  this  optimisation  problem  may  be  obtained  using  the  following  results.  Assume,  for  example, 
that  Cq  has  a  modified  PID  structure  of  the  form: 

C0  =  k0  +(kj  /  s)  +  (k2s  /  (1  +  sx))  (4.12) 

so  that  the  numerator  : 

C0/l  =  k0(l  +  sx)s  +  k1(l  +  sx)  +  k2s2  (4.13) 

and  the  denominator  : 

C0d=s(l  +  sx)  (4.14) 


Let  the  superscripts  r  and  i  denote  the  real  and  imaginary  parts  of  a  complex  function,  so  that 
C0n  =  Cgn  +  jC'0n  and  C0d  =  C’od  +  jC‘0d  .  The  controller  numerator  term  may  be  split  into  frequency 
dependent  components,  through  comparison  with  (4.13): 

2  2 

C0n  ( j®)  =  -k0(n  x  +  kt  -  k2u>  +  j( k0a>  +  kjox )  (4.15) 


and 


C()n( ~ 


x  +  kj-k2  co“  and  Cl0n  ( jm)  =  k0m  +  k,(ax 


(4.16) 


Similarly,  for  the  denominator  term: 


and  hence 


C0d(j(o)  =  -(o-‘z  +  j(o 


C'od  ( j ®  )  =  -® ' 2 1  and  C'0d  (j  m)  =  co 


(4.17) 

(4.18) 


If  the  solution  of  the  optimization  problem  is  to  be  found  by  iteration,  the  denominator  term  in  7}+  can 
be  assumed  to  be  known  and  the  minimisation  can  then  be  performed  on  the  numerator  (linear  terms).  Thus, 
to  set  up  this  problem  let, 

7}+  =  C0nLni  -  C0dLn2 


where 


Lni  ~  L]  /  (C0nL3  +  C0dL4)  and  Ln2  -  L2  /  ( C0nL3  +  C0dL4) 


Substituting  from  (4.1 1)  and  (4.19). 


t;  =  -  c‘Ai  -  «2 + C‘diln2  +  y(cUi + c'oAl  -  coA2  -  CM 


(4.19) 
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and  after  substitution  from  (4.16)  and  (4.18)  obtain 

K  {-0)2mn  1  - coCnX  +  j{coLrn]  - co2tL‘J)  +  kt  (r,  - (OTlinX  +  j(onLrn]  + 
+k2  [-or  Lni  -  jorl)n] )  +  orzl)n2  +  coFnl  +  j(a>~rFn2  -  o)Ln2  ) 


Tx+=\ 


The  real  and  imaginary  part  of  7)+  may  therefore  be  written  as:  7}4  =  7)1 '  +  /T/  '  and  it  follows  that, 


Write  a  vector  form  of  the  above  equations  as: 


where 


k0 

kj 

k2 


-L  =  Fx-L 


F(w)  = 


-coCco  xL'nl  +  L‘nI) 
m(L'n/-<mL,n,) 


(L'ni  -  c oxL'n] ) 

( (tixLrnl  +  Linl) 


2  Tr 

co  Lnl 

T‘ 

co  Lnl 


and  L( co)  - 


-co(  L'n2  +  co  xLrn2) 
co Lrn2-to2xlin2 


The  cost-function  can  be  optimised  directly  but  a  simple  iterative  solution  can  be  obtained  if  the  integral  is 
approximated  ( Yukitomo  et  al  1998  [23])  by  a  summation  with  a  sufficient  number  of  frequency  points 
{(0j,(02  ,...,(0N}  .  The  optimisation  can  then  be  performed  by  minimising  the  sum  of  squares  at  each  of  the 
frequency  points.  The  minimization  of  the  cost  term  Jq  is  therefore  required  where, 


N 

J0='L  (Fx- L)T (Fx- L)  =  (b- Ax)T (b- Ax) 

k=l 

where 


A  = 

~F(tOj)~ 

,  b  = 

L(a>i) 

,  x  = 

i 

_F(a>  N  )_ 

_L(  (oN)_ 

k2_ 

(4.20) 


(4.21) 


Assuming  the  matrix  AT A  is  non-singular  the  least  squares  optimal  solution  ( Noble  1969  [24])  follows  as: 


x  =  (ArAf1Arb 


(4.22) 


Lemma  4.1:  Restricted  Structure  LQG  Controller  Solution  Properties 

The  characteristic  polynomial  which  determines  stability  and  the  implied  equation  are  given  as: 

Pc  =  AC0d  +  BC0n 

AAqH0  +  BArG0  =  DcDf  (4.23) 
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The  minimum  value  of  the  cost-function,  with  controller  of  restricted  structure,  is  given  as: 

1 


Jnm  =  —  §  {TfTl  +TfTf  +®0}ds 


2xj 


2  71  j  D 


(Y0-F0Df)  ( Y0-F0Df ) 
- 1 - - - 

DcDcDfDf 


O0}t/5 


(4.24) 


Proof:  The  implied  equation  (4.23)  follows  from  (4.7)  and  (4.8)  by  multiplying  (4.7)  by  BAr  and  (4.8)  by 
AAq,  and  adding. 


4.1  Design  and  Robustness  Improvement 

The  LQG  controller  should  be  designed  in  such  a  way  that  it  is  consistent  with  the  restricted  controller 
structure  of  interest.  For  example,  Aq  should  approximate  a  differentiator  if  near  integral  action  is  required. 
In  fact  the  assumption  made  in  deriving  Theorem  3.1  was  that  the  controller  structure  is  compatible  with  the 
choice  of  error  weighting  and  if  1/Aq  includes  a  j  axis  pole  then  this  will  be  included  in  the  chosen  controller. 
In  fact  the  usual  situation  will  be  that  the  designer  decides  the  controller  should  include  integral  action  and  the 
weighting  ( 1/Aq )  will  be  chosen  as  an  integrator.  The  control  weighting  1/Ar  is  not  so  critical  but  if  for 
example,  a  PID  structure  is  to  be  used,  then  the  point  at  which  the  differential  (lead  term)  comes  in  can  help  to 
determine  the  Ar  weighting.  Clearly,  there  is  no  point  in  designing  an  LQG  controller  which  has  an  ideal 

response,  in  some  sense,  but  cannot  be  approximated  by  the  chosen  controller  structure.  Thus,  the  weightings 
should  be  selected  so  that  the  closed-loop  properties  are  satisfactory  but  taking  into  consideration  the 
limitations  of  the  controller  structure  required.  The  basic  concept  proposed  is  straightforward.  That  is,  in  the 
region  of  the  unity  gain  crossover  frequency  for  the  open  loop  system,  or  the  phase  margin  frequency, 
the  distance  \l  +  WCg\  should  normally  be  maximised.  This  requires  the  sensitivity  to  be  minimized, 
particularly  in  this  sensitive  region.  By  costing  the  sensitivity  directly  a  mechanism  is  provided  to  improve 
robustness  (Horowitz  1979  [18])  but  there  are  some  subtleties  to  address: 

1.  The  weighting  Pc  needs  to  be  increased  from  zero  where  performance  is  presumably  maximized 
(the  LQG  cost  is  optimized)  up  to  a  level  where  robustness  is  adequate  and  performance  still  acceptable. 

2.  If  pure  sensitivity  costing  is  required  Pc  could  cancel  the  combined  noise  dynamics  ®  ^  =  Yj  Yj  which  is 

unrealistic.  The  alternative  is  to  make  the  model  W  =  A~1Cp  large,  relative  to  the  other  noise  terms  and 

also  a  constant  ( C  =  pA  say)  and  this  will  introduce  a  fictitious  stochastic  term  affecting  the  noise  and 

disturbance  rejection,  and  reference  tracking,  properties.  The  size  of  the  scalar  p  also  therefore  involves  a 
compromise. 

3.  The  weighting  and  frequency  shaping  effects  are  only  important  in  the  decade  above  and  a  little  below  the 
crossover  frequency  referred  to.  The  shaping  might  therefore  be  introduced  by  a  weighting  that 
approximates  an  ideal  window  function  but  this  increases  the  order  of  the  weighting  term. 

The  above  design  choices  and  trade-offs  detract  from  the  approach  but  the  prize  is  quite  important  and  worthy 
of  the  effort.  That  is,  the  provision  of  a  tuning  variable,  or  variables,  in  a  cost  index  where  the  robustness 
properties  of  an  optimal  controller  can  be  manipulated  and  traded  against  performance/stochastic  properties. 
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5.0  COST-FUNCTION  EVALUATION  AND  CONTROLLER  BENCHMARK  COST 
COMPUTATION 

The  minimisation  of  the  cost-function  can  be  shown  to  be  equivalent  to  the  minimisation  of  the  variance  of  the 
adjoint  output  signal  (see  Figure  2). 

=  D:-\AlRcXu-B;QcnAq\ )  (5.1) 


A 


Dc  1  BjQcnAq1 


'Dc^lRcnAr1 

I _ 

▲ 


d 


n 


Figure  2:  Block  Diagram  representation  of  Benchmarking  Cost. 


The  polynomial  operator  version  of  the  control  spectral  factor  may  be  defined  using: 

7*7  =W*QW  +  R  =  A*\B*Q  B,  +  A*R  A,)A~l 

c  c  c  c  v  1  z-scn  1  1  cn  1 '  c 

where  Dc  is  a  polynomial  operator  that  satisfies: 

D*D  =  B*Q  B,+A:R  A, 

c  c  1  zZ'Cn  1  1  cn  1 

and  Yc  may  be  written  as:  Yc  -  DCA~1  where  Aj  =  AAg  and  Bj  =  BAr,  and  Ac  =  AAgAr  Assume  the  existence 
of  the  solutions  (G0,Fq),  (H0,Fq)  of  the  equations  (4.7)  and  (4.8).  These  equations  can  be  added,  after  first 

multiplying  by  Dj^Bj  and  DJ1  A} ,  respectively,  to  obtain: 
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Dc(G0Df1B1+H0D-f1A1)  +  F0(AAqD-f1B1  -BArDf' A,) 

=  (Z;QcA+4RcA)  (5.2) 

but  the  last  two  of  the  terms  in  (5.2)  add  to  zero.  The  implied  diophantine  equation  then  follows  from  (5.2)  by 

* 

dividing  by  Dc  : 

G0Bj  +  HqAj  =  DcDf  (5.3) 


The  variance  of  such  a  signal  (f)  =  Dc  1  (  A* R cn  A r '«  -  B*QcnAq  'e0)  can  be  shown  (using  ParsevaTs  theorem) 
to  be  the  same  as  the  signal  ( Grimble  1984  [18]): 

4  =  D;lDl+  =  D:\AXXu  - BlQcnA-q\) 


The  solution  of  the  LQG  optimal  control  problem  therefore  involves  the  minimization  of  the  signal  shown 
(j)cj  in  Fig.  2.  This  is  valuable  in  providing  a  link  to  the  benchmarking  techniques  for  minimum  variance 
controllers  that  are  well  established.  The  diophantine  equation  (4.7)  may  be  written  in  the  form: 


Dc  H 


Fn 


+  -RT  =  D~lB*  QcnAq\AlDf 


D 


Recalling  that  /  =  A  1  Dj  s  the  equation  above  may  be  used  to  obtain: 

A~\f 

D  AA„  D  c  q 


From  this  last  result  one  interpretation  of  the  term  F0D~]  e ,  which  has  the  same  variance  as  part  of  the  cost 
term  Jy,  (defined  in  (5.6))  is  one  component  of  /  in  the  signal  ^  in  Fig.  2.  In  minimum  variance 
benchmarking  problems  this  term  is  of  the  simpler  form  Fq  e  and  the  coefficients  of  Fq  are  available  from 
simple  test  results  involving  the  signal / (considered  in  the  discrete  time  case  by  Huang  and  Shah  1999  [30]). 

The  value  of  the  cost-function  and  of  the  error  and  control  signal  variances  can  be  evaluated  using  the  solution 
of  the  equations  in  Theorem  3.1.  The  following  theorem  provides  the  expressions  for  the  cost  values,  for  both 
the  case  of  any  (suboptimal)  controller  and  for  the  optimal  solution. 

Theorem  3.2: 

The  performance  criterion  (3.1)  may  be  evaluated  for  any  linear  controller  C0  =  Cq^Cqh  as: 

J  =  J0  +  Ji  +  J2  (5.4) 

where 

Jo  iT^}ds  (5-5) 

2*Jd 
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•/,  = 


2  n  j  \ 


(  Y  -  F  D*  \ 

-*0  1  0  ^  f 

f  Yo  ~  F*D)  \ 

D  *,D  * 

V  f  c  J 

d  ;  d  * 

V  f  c  J 

yds 


(5.6) 


Jl  2  nj 


<D  ,  <D  , 

lp  fp  V ds 


c  ff  z-' c  nn  y  V*V  V* 
D  {  1  c1  c  1  f1  f 


where 


w  = 


\(H()All()Aq  +  X0ArB)C0n  ( G0A/l0  X0A)ArC0d\ 
AiO^gA-i^Od  +^Q/i) 


(5.7) 


(5.8) 


The  variances  of  the  error  and  control  signals  may  also  be  computed,  for  any  controller,  as: 


J.  =  A~A  U7«>  +  A-f  {«„„  -os +«>„)}<& 

2 KJD  2nj  d 


(5.9) 


where 


*4  =  t — :§  {TuT*}ds 
2k  J  D 


Te  ~  C0d  ( AC0d  +  BC0n  )  Df  and  Tu=  C0Te 


(5.10) 


(5.11) 


The  optimal  value  of  the  cost  term  Jq,  for  the  full  order  optimal  controller,  is  Jf)mjn  =  0  .  The  minimum-cost 
in  the  full-order  optimum  case  can  be  obtained  in  the  alternative  forms: 

A  min  —  J 1  +  J 2  (5-12) 


where 


and 


Anin  d q  min  T  A- min  dp  min 


Ti.il.  =  A-f  {smJ,QX,sL +®.a-(s.  ,n  +0&®,,,)}* 

2-K  J  D 

Jr™=A-§  =  A- f  {T-C,  )* 

2 kjd  2nj  d 

T«... = A-nf  {s«*Y,pX,sL)&  =Amf  { 

2 njDy  2n ] d  y 


_B,.(G0An0-X0A)  _Bp(H0An0Aq-X0ArB) 

d-  v  min  —  5  1  p  mm 


DCU o 


dqDcU0 


(5.13) 

(5.14) 

(5.15) 

(5.16) 

(5.17) 
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The  optimal  values  of  the  variances  for  the  error  and  control  signals  for  the  full-order  optimal  controller 
follow  as: 


where 


and 


J. 


=  —'—(6  \t  T*  .  +0  -( s  +5*.  )®  Ids 

j  J  (_  emm  ermn  nn  v  rrun  min/  nn  J 


}ds  =  (L,nC,n)^ 

D  An  j  D 


V  .  — 

^min 


Ao(H0An0  A  +  X0ArB ) 


M  ■  = 

min 


A0Ar(G0An0-X0A) 

DcDf 


_A,.(G0An0-X0A)  _H0An0Ag-X0ArB 

1  u  min  T  T  5  1  e  min 


DCU  o 


(5.18) 

(5.19) 

(5.20) 

(5.21) 

(5.22) 


The  existing  system  might  have  a  classically  designed  controller,  of  more  conventional  structure  (like  a  P1D 
controller),  and  although  it  is  of  some  value  to  compare  this  cost  J  for  this  classical  controller  with  the  best 
that  can  be  achieved  ( Jmin  =  Jj  +  J^),  this  is  often  an  unfair  benchmark  comparison.  The  reason  is  of  course 
that  the  restriction  on  the  controller  structure  and  the  order,  may  make  the  absolute  minimum  Jm[n  way  below 

what  is  achievable  using  the  best  tuned  classical  controller.  A  more  appropriate  benchmark  figure  is 
therefore  to  compare  the  actual  cost  with  that  of  the  optimal  cost  assuming  the  controller  structure  is  fixed. 
If  the  controller  has  a  specified  limited  structure,  the  minimum  of  the  cost  term  Joinin  be  non-zero. 
The  minimum  value  of  the  criterion  (3.1)  for  this  restricted  controller  will  be  obtained  as 
J  restrict  =  ^Omin  +  J  min  ■  b°r  an  unconstrained  solution  the  minimum  is  achieved  when  7j+  =  0  and  the 
minimum  of  Jq  (denoted  Jomin )  is  zero.  The  increase  in  cost  which  occurs  by  restricting  the  controller 
structure,  can  be  obtained  as, 

=  (TTV  (5.23) 

D 


The  Controller  Performance  Index  will  be  defined  as  the  ratio  of  the  minimum  possible  value  of  the  cost 
function  (3.1)  to  the  actual  value  and  the  CPI  lies  between  0<k<1  .  If  K  is  close  to  unity  the  system 
provides  little  opportunity  for  improvement.  If  the  CPI  is  close  to  zero  retuning  is  recommended.  This  scalar  is 
similar  but  not  the  same,  as  the  assessment  measure  introduced  by  Desborough  and  Harris  (1992  [29]) 

^  ~  ^min  /  (Jmin  Jq  )  ~  ^  ~  Jq  /  (Jmin  Jq  )  (5-24) 


RTO-EN-SCI-142 


2-17 


Design  and  Performance  Assessment  of 

Robust  Restricted  Structure  Optimal  Control  Systems 


ORGANIZATION 


6.0  ROBUST  CONTROL  DESIGN  EXAMPLE 

To  illustrate  the  effect  of  the  robustness  weighting  element  and  the  fictitious  robustness  signal  {p(t)}  a  simple 
example  is  considered.  The  system  models  may  be  listed  as: 


Coprime  system  model: 


W=^-  =  - 


1000(s  +  2)(s  +  6) 


A0  (s  +  0.7)(s  +  3.9)(s  + 100)(s2  +  2  £y0s  +  co 02) 
where  c  =  0. 1  and  co0  =  10.  Also  write 


Disturbance  model: 

where 

Reference  model: 

where 

Noise  model: 

where 


w-l-ML o 

A  A0U0 


where  Uq  =  (s+3.2)s 


1000 


A  (s  +  100)(s2  +2^co0s  +  o>q  )(s  +  3.2)5 
Q  =1000(5  +  0.7)(s  +  3.9) 


IF. -4- 


1 


A  ( s 2  +  2  %cl>qS  +  col  )s 

E  =  (s  +  0.7)(s  +  3.9)0  +  3-2)0  +  100) 


=  -  01 


"  A  O  +  100) 

Cn  =  0. 10  +  0.7)0  +  3.9)02  +  2  ga>0s  +  ca\  )0  +  3.2  )s 


Fictitious  robustness  signal: 


W„ 


P 

O  +  ioo) 


Cost  Function  Weightings 

The  cost  function  weightings  may  be  defined  as: 

_  (0.01s +  l)(-0.0  Is +  1) 
e  ~  (s  +  10“6)(-s  +  10"6) 

Rc  =  (10s  +  l)(-10s  + 1) 

„  (0.01s +  1X-0.0  Is +  1) 

Pc  =  Pi - 7 - 7— 

(s  +  10“6)(-s  +  10“6) 
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Results 

The  frequency  responses  of  the  different  system  models  are  shown  in  Fig.  3.  The  system  is  low  pass  with  a 
resonant  subsystem  and  the  disturbance  model  includes  an  integrator.  The  measurement  noise  model  only  rolls 
off  at  high  frequencies.  Consider  first  the  full  order  optimal  case  and  the  use  of  a  large  p  — 1 000  ,  then  as 

px  varies  the  unit  step  responses  of  the  closed  loop  system  are  as  shown  in  Fig.  4.  This  represents  the  case 
where  there  is  a  large  fictitious  disturbance  model  Wp  but  where  the  robustness  weighting  p  varies  between 
0  to  1000.  The  faster  responses  occur  as  p  increases,  since  the  effect  is  related  to  that  when  Qc  increases. 
The  corresponding  closed-loop  frequency  responses  are  shown  in  Fig.  4. 


Figure  3:  Bode  Frequency  Responses  for  the  System,  Disturbance  and  Noise. 


Figure  4:  Closed  Loop  System  Unit  Step  Responses  for  Fixed  p  =  1000  and  Varying  Weighting  p^. 
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Figure  5:  Bode  Diagram  of  Closed  Loop  Responses  for  p  =  1000  and  Changing  Value  of  p-j. 
Restricted  Structure 

The  unit  step  responses  are  compared  in  Fig.  6  for  the  full-order  and  restricted  controller  designs.  The  case 
p-  0  and  px- 10  was  considered  and  the  corresponding  controller  and  sensitivity-function  frequency 
responses  are  shown  in  Figs  7  and  8.  The  restricted  structure  control  is  particularly  good  from  an  overshoot 
perspective.  However,  one  reason  is  the  higher  controller  gains  at  high  frequencies  for  the  restricted  structure 
control  law.  The  computed  controllers  were  obtained  as: 

Optimal  Full-Order  Controller: 

c  _  0.493426(5  +  0.2205695)(s  +  0.7)(s  +  3.199631)(s  +  3.9)Q2  +  1.999939s  +  99.99969 )(s  +  100) 

°  'S  (5  +  0.00000 1)0  + 1.09342 X  10“6)(s  +  1.008371)(s  +  3.866217)(s2  + 14.817935  +  65.05414) 

x(52  +7.5208615  +  166.0202) 


Optimal  Restricted  Structure: 


C0(s)  = 


-0.064124852  -  0.358762b  +  1.208584 
5(0.25  +  1) 
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Figure  6:  Comparison  of  Closed-Loop  Unit  Step  Responses  of  Full  and 
Restricted  Structure  Control  Designs. 


Figure  7:  Bode  Comparison  of  Controller  Frequency  Responses. 
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Figure  8:  Bode  Comparison  of  Sensitivity  Function  Frequency  Responses. 

The  step  responses  shown  in  Fig.  9  are  for  the  case  p-  P\  =100.  The  results  are  much  faster  and  the 
restricted  structure  design  is  again  good,  relative  to  the  full-order  solution. 


Figure  9:  Comparison  of  Unit  Step  Responses  for  Full-Order  and 
Restricted  Structure  Control  Designs. 
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Example  Conclusions 

The  example  reveals  that  the  robustness  weighting  terms  and  the  fictitious  robust  costing  signal  {p(t)} 
certainly  affects  the  overshoots  which  represent  a  measure  of  robustness,  both  on  closed-loop  frequency  and 
time  responses.  The  tuning  variables  p  and  pl  affect  the  robustness  of  this  minimum-phase  open  loop  stable 
system  in  much  the  way  expected.  However,  the  alteration  of  robustness  properties  is  not  a  straightforward 
matter,  since  any  values  of  p  and  /?,  above  zero  will  cause  a  measure  of  sub-optimality  in  stochastic 

(LQG  cost)  terms.  The  most  surprising  results  were  the  very  good  results  obtained  for  the  restricted  structure 
control  designs.  The  explanation  was  the  higher  high  frequency  gains  employed  that  reduced  the  peaks  on  the 
sensitivity  function  frequency  responses.  In  this  problem  changes  in  the  measurement  noise  model  again  did 
not  have  a  large  effect.  However,  results  that  were  not  shown  were  obtained  for  a  coloured  measurement  noise 
model  with  a  peak  in  the  low  frequency  range.  In  this  case  the  controller  gains  are  significantly  reduced  and 
this  slows  the  speed  of  response  of  the  system  and  the  overshoot  increases  markedly. 


7.0  CONCLUSIONS 

The  robustness  of  full-order  and  restricted  structure  optimal  control  problems  was  considered  for  continuous¬ 
time  linear  systems.  The  emphasis  was  on  the  improvement  of  robustness  by  adding  a  sensitivity  costing  term 
in  the  cost  index  and  by  introducing  a  fictitious  disturbance  model  Wp.  The  effect  of  measurement  noise  was 

also  investigated  and  this  was  introduced  in  the  feedback  system  model.  The  robustness  weighting  acts 
directly  on  the  sensitivity  function  and  may  improve  robustness  margins  but  other  properties  will  probably 
deteriorate  like  the  measurement  noise  rejection  properties.  If  robustness  is  more  important  than  stochastic 
properties  then  attention  would  turn  to  Hx  cost  minimization  and  many  of  the  ideas  presented  above  would 
apply  ( Grimble ,  1986  [20]).  However,  such  an  approach  is  readily  embedded  in  the  usual  mixed  sensitivity 
H ^  design  problem. 
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